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ABSTRACT 

We present the high resolution spectra of the D type symbiotic stars V1016 Cygni and HM Sagittae 
obtained with the Bohyunsan Optical Echelle Spectrograph (BOES), and investigate the double- 
peaked asymmetric profiles of the Raman scattered O VI 6825. This feature is formed through 
Raman scattering of O VI 1032 by atomic hydrogen with small scattering cross section of ~ 10 -23 cm 2 , 
requiring a specific condition of coexistence of a highly ionized emission nebula and a thick neutral 
region. By adopting a wind accretion disk model, we assume that the O VI emission region is described 
by a Keplerian thin disk. The Raman scattering occurs in a neutral region near the giant, taking in 
the form of a slow stellar wind, part of which is ionized by the strong UV radiation from the hot 
white dwarf. Using a Monte Carlo technique, we compute the line profiles that are modulated by 
the slow spherical stellar wind from the giant component with the ionization front approximated by a 
hyperboloid. In order to account for the asymmetry and the existence of a central dip in the profiles, 
we add an O VI resonance scattering region between the hot white dwarf and the giant star which 
hinders the incidence of slightly blue O VI photons upon the H I region. Overall good fits to the 
observed data are obtained from our model, which lends support to the accretion disk emission model 
in these objects. The best fitting parameters for V1016 Cyg are v = 30 km s _1 , Voo = 11 km s _1 , 
and v c — 10 km s _1 , where v , Voo and v c are the velocity of the outer disk rim, the terminal velocity 
of the giant wind, and the velocity component of the resonance scattering O VI region along the binary 



axis, respectively. Similar fitting parameters v a — 27 km s 



10 km s and v c 



9 km s 



are obtained for HM Sge. We also investigate the effect of a hot spot in a disk that is well known 
in accretion disks in cataclysmic variables. However, the introduction of a hot spot in our Keplerian 
disk model failed to improve the overall profile fitting quality significantly. Brief discussions about 
our profile analysis in relation to bipolar morphology and accretion processes are presented. 
Subject headings: binaries : symbiotic — accretion, accretion disks — line : profiles — scattering — 
stars (individual V1016 Cyg, HM Sge) 



1. INTRODUCTION 

Symbiotic stars are spectroscopically characterized by 
prominent emission lines with TiO absorption bands that 
are typical of a giant. They are generally known to be 
wide binary systems consisting of a hot white dwarf and 
a giant star (e.g. Kenyon 1986). They are usually di- 
vided into 'S' and 'D' type symbiotic stars, where 'D' 
type systems exhibit infrared excess indicative of warm 
dust with temperature T ~ 10 3 K. Many 'D' type sym- 
biotic stars are believed to contain a Mira type giant as 
cool component. Their orbital separation of the hot and 
cool components may lie in the range of 10-100 AU, from 
which we expect the orbital periods of order 10 2 years 
(Whitelock 1987). 

The giant component may lose mass in the form of 
slow stellar wind, a significant fraction of which can be 
photoionized by the hot white dwarf component. Taylor 
& Seaquist (1984) investigated the ionization structures 
in order to interpret their radio observations of symbi- 
otic stars. Many prominent emission lines in symbiotic 
stars can be excellent diagnostic tools to study the phys- 
ical conditions around the white dwarf component (e.g. 
Murset et al. 1991). 

It has been proposed by a number of researchers that 
some fraction of wind material from the giant compo- 
nent may be gravitationally captured by the hot compo- 



nent (e.g. Sokoloski Bildsten & Ho 2001). The mass 
loss rate from the giant component may range from 
10~ 4 — 10~ 7 M Q yr _1 depending on the stellar evolution 
stage (e.g. Iben & Tutukov 1996, Luthardt 1992). Wind 
accretion in binary star systems was investigated using 
the smoothed particle hydrodynamics (SPH) by Theuns 
and Jorrissen (1993), who concluded that an accretion 
dis can be formed in their isothermal flow models. Mas- 
trodemos & Morris (1998) also adopted an SPH method 
to show that an accretion disk can be formed around the 
white dwarf component. There are only limited studies 
about the existence of accretion disks and the insuffi- 
cient numerical resolution prevents one from deducing 
the physical properties of an accretion disk in symbiotic 
stars. Symbiotic stars being considered as a candidate for 
type la supernova progenitors, the mass transfer process 
in symbiotic stars can also be a very important subject 
in other branches of astronomy. 

About a half of symbiotic stars are known to exhibit 
very broad and mysterious emission features around 6825 
A and 7088 A that were suspected to originate from 
highly ionized species (Allen 1980). These features were 
finally identified by Schmid (1989), who proposed the 
Raman scattering origin. According to him, the Raman 
scattering process starts with an O VI 1032 line pho- 
ton incident upon a hydrogen atom in the ground state, 
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which can be converted into an optical 6825 photon as 
a result of an inelastic scattering off the hydrogen atom 
that is left in its excited 2s state. From an incident O VI 
1038 line photon, we obtain a Raman photon at around 
7088 A. The Raman scattering nature of these features 
is strongly supported by many observational studies in- 
cluding spectropolarimetric observations by Harries & 
Howarth (1996) and contemporaneous UV and optical 
observations by Birricl et al. (2000). 

The wavelength A Q of the inelastically scattered radia- 
tion is related with the incident wavelength Aj and A L ya 
of hydrogen Lya by 

A c 1 = \ 1 — \ a 1 . (1) 

Differentiation of this equation leads to 



In the case of Raman scattering involving O VI 1032, 
A /Ai = 6.6. This provides two important characteristics 
in the profiles of Raman scattered O VI 6825 that they 
are broader than their parent emission lines by a factor 
of X /Xi and that they are mainly determined by the 
relative motion between the neutral scatterers and line 
photon emitters. Therefore, the profiles are little affected 
by the direction of the observer's line of sight. This is 
due to the inelasticity of Raman scattering and noted to 
be one of the most important spectroscopic properties 
retained by Raman scattered features (e.g. Nussbaumer, 
Schmid & Vogcl 1989, Schmid 1989). In this respect, 
Raman scattered lines provide a very interesting oppor- 
tunity to study the accretion process involving a slow 
stellar wind, which is quite unique to symbiotic stars. 

Thus far, Raman scattered lines are found for O VI 
1032, 1038 resonance doublets and He II emission lines 
near H I Lyman series in many symbiotic stars and in a 
few young planetary nebulae (e.g. van Groningcn 1993, 
Pcquignot et al. 1997, Groves et al. 2002, Lee et al. 
2003, Birriel 2004, Zhang et al. 2005, Lee et al. 2006b). 
Raman scattered O VI 6825 from O VI 1032 is particu- 
larly strong in a number of symbiotic stars, for which a 
refined profile analysis is feasible. Furthermore, Raman 
scattered features consist of purely scattered photons, 
not mixed with the direct emission, which makes them 
characterized by strong linear polarization. 

Fairly intensive spectropolarimetric studies of Raman 
scattered O VI lines were carried out by Harries & 
Howarth (1996). They found that the Raman scattered 
O VI 6825 and 7088 show double or triple peak profiles 
and that they are strongly polarized with complicated 
structures. Schmid (1996) used a Monte Carlo method to 
calculate the profiles and polarization of these scattered 
features formed in a slowly expanding neutral region that 
mimics the stellar wind around the giant component. A 
similar study with a slightly different velocity law of the 
slow stellar wind was performed by Lee & Lee (1997b). 

In these studies, no kinematic structure is assumed 
for the emission region, which is treated effectively as 
point-like without any complicated internal structures. 
Instead, the stronger red part of Raman scattered O VI 
6825 is mainly attributed to scattering at the receding 
part of the giant stellar wind that has a bigger geometric 
covering of the O VI emission region than the approach- 
ing part does. However, in these models, one has to 



invoke a rather high velocity scale > 50 km s _1 of the 
stellar wind from the giant in order to explain the ob- 
served Raman profiles. 

Lee & Park (1999) adopted an accretion disk emis- 
sion model and proposed the origin of the multiple peak 
profiles with strong linear polarization. In their picture, 
the accretion flow can be affected significantly by the 
binary orbital motion, which results in an asymmetric 
matter distribution around the white dwarf causing the 
stronger red part of the Raman features. However, their 
study lacked detailed calculations that can be directly 
compared with the observed profiles. 

Both V1016 Cyg and HM Sge are classified as symbi- 
otic novae or slow novae. V1016 Cyg underwent a slow 
nova-like outburst in 1965, which may be interpreted to 
be a thermonuclear runaway on the surface of the white 
dwarf (Fitzgerald et al. 1966). A similar nova-like op- 
tical outburst occurred in HM Sge in 1975 (Dokuchacva 
1976). These two symbiotic stars are known to exhibit 
fairly strong Raman scattered O VI 6825 with character- 
istic double peak profiles. In particular, the appearance 
of Raman scattered 6825 in HM Sge was reported by 
Schmid et al. (2000) in their spectra obtained with the 
William Herschel Telescope (WHT) in 1998. 

Both V1016 Cyg and HM Sge are D-type symbiotic 
stars, containing a Mira type giant as cool component. 
Thus far, no detailed information is available about the 
orbital parameters of V1016 Cyg and HM Sge. Schild 
& Schmid (1996) measured the rotation of position an- 
gle in the polarized Raman 6825 with an estimation of 
the orbital period of ~ 80 yr in V1016 Cyg. However, 
Schmid (1998) noted that this estimate may be severely 
affected by the change of the nebular structure in this 
system and the orbital period may exceed 100 years (e.g. 
Schmid 2001). Using the Hubble Space Telescope (HST) 
data, Brocksopp et al. (2002) measured the angular sepa- 
ration of 42.4 mas between the white dwarf and the Mira 
of V1016 Cyg, proposing the physical binary separation 
of <~ 84 AU with their adopted distance of 2 kpc. Ac- 
cording to the HST study of HM Sge by Eyres et al. 
(2001), the separation of the Mira and the white dwarf 
is <~ 50 AU with the adopted distance of ~ 1.2 kpc. 

In this paper, we present our spectra of these Raman 
features in the symbiotic stars V1016 Cyg and HM Sge 
secured with the Bohyunsan Optical Echelle Spectro- 
graph (BOES), and make quantitative comparisons with 
the profiles computed by a Monte Carlo technique adopt- 
ing an accretion disk emission model. 

2. OBSERVATION AND DATA 

We observed V1016 Cyg and HM Sge on the nights of 
2005 November 7 through 9 with the Bohyunsan Optical 
Echelle Spectrograph installed on the 1.8 m telescope at 
Mt. Bohyun. The detector was a 2kx4k E2V CCD with 
pixel size 15 /im x 15 /im, in which about 80 spectral 
orders covering ~ 3700 A to ~ 10,000 A are recorded 
in a single exposure. The 300 fx optical fiber was used 
to yield the spectroscopic resolving power of R = 30, 000 
with an aperture field of 6.4 arcseconds. A more detailed 
description of BOES can be found in Lee et al. (2006b) 
and in Kim et al. (2002). The exposure times are from 
7000 s and 3600 s for V1016 Cyg and HM Sge, respec- 
tively. The data have been reduced following standard 
procedures using the IRAF packages. 
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Parts of our spectrum around 6825 A of V1016 Cyg 
and HM Sge are shown in Fig. 1. In the figure, the 
line flux is normalized by the local continuum in order 
to focus on the profile shape in the current work. The 
two profiles are overall similar to each other. The 6825 
feature of HM Sge is much weaker and therefore of poorer 
quality than that of V1016 Cyg, which is indicated by the 
larger scatter in the data. In both systems, the red part 
is relatively stronger than the blue part, which is noted 
by many previous researchers (e.g. Harries & Howarth 
1996, Schmid et al. 1999). 

There exists a narrow emission line at around 6820 A in 
HM Sge. The presence of this narrow emission feature in 
HM Sge was also noted in the spectropolarimetric stud- 
ies by Schmid et al. (2000). They concluded that this 
emission line is not related with Raman scattering on 
the ground that it is unpolarized and narrow. We ignore 
this feature in the subsequent analysis, by subtracting it 
by a single Gaussian /(A) = /oexp[(A — A C ) 2 /AA 2 ] with 
/o = 0.43, A c = 6820.09 A and AA = 0.46 A. 

In the case of V1016 Cyg, the local maxima are found 
at 6823.5 A and 6830.8 A. The separation of these two 
peaks corresponds to 48 km s _1 in the velocity space of 
O VI 1032 emitters. The FWZI (full width at zero inten- 
sities) is 166 km s -1 . It is noted that there is an asym- 
metrically extended red wing that may end at around 
6842.7 A. 

The peaks for HM Sge are located at 6823.3 A and 
6829.4 A, for which the separation in velocity space of 
O VI is 41 km s _1 . This is significantly smaller than 
that for V1016 Cyg. The FWZI is 161 km s -1 , which 
is comparable to that of V1016 Cyg. In contrast with 
V1016 Cyg, it appears that there is no (or significantly 
weak) extended red wing in HM Sge. However, this 
should be taken with the caution that the data quality 
of HM Sge is much poorer than V1016 Cyg. 

The relative kinematics between the O VI emission re- 
gion and neutral H I scattering region in both V1016 Cyg 
and HM Sge is characterized by velocity scales ranging 
from 20 km s" 1 to 80 km s _1 . If the FWZI of the pro- 
files is attributed to the kinematics of the giant wind as 
in the models of Schmid (1996) and Lee & Lee (1997b), 
then the wind terminal speed of 80 km s _1 should be re- 
quired, which is much faster than the escape velocity of 
a giant star. 

3. MONTE CARLO LINE PROFILE ANALYSIS 
3.1. Accretion Disk Model 

In this work, we assume that the wind accretion occurs 
in both V1016 Cyg and HM Sge to form an accretion 
disk around each white dwarf component. According to 
Model 3 of the SPH studies by Mastrodemos & Morris 
(1998), an accretion disk with size ranging 0.4 - 0.8 AU 
is formed in a binary system containing a white dwarf, 
where the binary orbital period is 18 years. We consider 
that the O VI emission region is limited to a part of the 
whole accretion flow, which is approximated to be Kep- 
lerian. For simplicity, we assume that the O VI emission 
region is confined to an annular region characterized by 
the inner radius Ri and the outer radius R Q . 

In this work we assume that each white dwarf com- 
ponent has a mass of 0.7 M©. Noting that the Raman 
scattered O VI 6825 in both V1016 Cyg and HM Sge 
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Fig. 1. — Parts of spectra around 6825 A of the symbiotic stars 
V1016 Cyg and HM Sge obtained with the Bohyunsan Optical 
Echelle Spectrograph (BOES). The vertical scale is normalized by 
the local continuum. The peak separation is ~ 48 km s" 1 for 
V1016 Cyg and ~ 41 km s" 1 for HM Sge. The full widths at 
zero intensities are 166 km s — 1 and 161 km s — 1 for V1016 Cyg 
and HM Sge, respectively. 

extends about Av ~ 160 km s _1 , we infer that the O VI 
emission region extends from ~ 0.05 AU to less than 1 
AU. 

With these ingredients we prepare a characteristic dou- 
ble peak profile of O VI 1032, we first consider the disk 
component with the emissivity e(r) that depends only on 
r through a power law 



e(r) = Ar a , Ri < r < R Q 



(3) 



where A is a constant. This functional dependency of 
emissivity is achieved in our Monte Carlo code by speci- 
fying the location r of the initial photon from a uniform 
random number r p in the interval (0, 1) with the pre- 
scription 



Ru 



Ri 
Ru 



a+2 




l/(a+2) 



(4) 



With another random number in the interval (0, 2ir), 
we locate the point of generation of the incident photon 
P(r,<f> = Then the photon is assigned a Doppler 

factor DF along the propagation direction 



DF = ki • v(P)/c 



(5) 



where v(P) = J GM ™ D (- sin (f>x + cos^y) is the local 

Keplerian velocity at P. 

According to Eq. ([2]), the observed wavelength of Ra- 
man scattered O VI is determined by the relative motion 
between the O VI emitter and the H I scatterer and quite 
insensitive to the observer's line of sight. Furthermore, 
the orbital parameters of V1016 Cyg and HM Sge are 
only poorly known. In view of these facts, in our Monte 
Carlo code, we collect all the photons scattered inelas- 
tically off hydrogen atoms around the giant component. 
However, a more refined study should include the effect 
of the observer's line of sight with full care. 

3.2. Velocity Modulations from a Slow Stellar Wind 
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The direct O VI 1032 line photons are scattered in a 
neutral scattering region, which may be expanding with 
a velocity ~ 10 km s _1 . Schmid (1996) investigated the 
line profile formation of Raman scattered O VI in an 
expanding slow stellar wind using a Monte Carlo tech- 
nique (see also Harries & Howarth 1997). A similar study 
was also conducted by Lee & Lee (1997b) using a Monte 
Carlo method with a slight different kinematical velocity 
law for the slow stellar wind. 

Adopting also a Monte Carlo technique, we extend the 
work of Lee & Lee (1997b) to incorporate the accretion 
disk emission model considered in the previous section. 
The velocity of the slow and spherical stellar wind v#(r) 
at the position r is given by 



,(l-R*/rY 



(6) 



where is the terminal wind velocity and the origin 
of the coordinate system coincides with the center of the 
giant. For simplicity of the current investigation, we set 
the power j3 = 1 of the velocity law. We introduce a new 
dimcnsionlcss vector defined by 



p = r/R*. 



(7) 



Assuming a constant mass loss rate M, the density nr) 
is described by 

n(r) =n Q p- 2 (l-p- 1 )-\ (8) 

where no — M / '4ir pR 2 m p v ^ with m p and p being the 
mass of a proton and the mean molecular weight, respec- 
tively. 

In the Monte Carlo calculation, it is essential to locate 
the next scattering position p2 for a given starting po- 
sition p\ in the direction k^, the unit wavevector of the 
photon under consideration. In order to accomplish this, 
an optical depth t\2 is generated using a uniform random 
number r tau between and 1 by 



T%2 



In r ta 



(9) 



In terms of the physical distance s\i corresponding to t\i 
from the starting position p\ along the direction k;, the 
new scattering position is expressed as 

Pi = Pi + si 2 kY (10) 
The new scattering position p\ is related with m by 



T"l2 



n(r)a to t dl, 



(11) 



pi 



where a to t is the sum of the cross sections for Rayleigh 
and Raman scattering. Lee & Lee (1997a) proposed the 
value atot — 42ctt — 2.8 x 10 -23 cm 2 with ar = 6.6 x 
10~ 25 cm 2 being the Thomson scattering cross section, 
which we adopt in this work. 

We introduce the representative optical depth tq de- 
fined by 



t q = n R*a tot = 8.8M_ 6 i?*i 3 v, 



00IO7 



(12) 



where M_ 6 = M/(10" 6 M Q yr" 1 ), R* 13 
i?»/(10 13 cm) and v^w = u oo /(10 km s" 1 ). The 
exact choice of tq and a tot affects the ionization 
structure and hence the exact profile shape. However, 
the ionization structure is also dependent on the UV 
radiation from the hot component. In view of this, we 



set tq — 2.5 in our simulations, in order to reduce the 
parameter space of investigation. 

We also define b as the impact parameter of the photon 
path with respect to the giant center, and the parameter 
s measures the physical distance along the photon path. 
For the sake of convenience, we set s — at the foot of 
the perpendicular from the giant center to the photon 
path, so that 

s = ±^/p 2 ~ V, (13) 

where negative s is obtained when k^ • p\ < 0. In terms 
of these parameters, the integral (fTTj) can be rewritten 



T12 



n(s) a to t ds = T Q I 

Si J pi 



dp 



In terms of a new function u(p) defined by 
1 



(14) 



u(p) 



jp-i + vv^n (15) 



the integral in Eq. (fl4| can be explicitly expressed as 
2r 



T12 = 



Using the inverse relation of Eq. (fT5|) 



[tan u(pi) — tan u(pi)]. (16) 



(17) 



p = 6coshln[(\/6 2 - 1m + l)/6] 

this relation is directly inverted to obtain P2 in terms of 
T12 and p\. 

A more complicated situation is obtained when b > 
l,kj • ri < 0, in which case we consider 77,, the optical 
depth to the point of impact, defined by 



n = 



2r 



If T12 > tj, we consider 
ri2 = n + 
from which we obtain 

u(p2) = tan 



tan 1 u\ — tan 1 



6-1 



1/2- 



n{p) a tot ds, 



\/b 2 - 1(7-12 - n) 
2m 



(18) 



(19) 



(20) 



With p2, S12 is determined, and the location of the next 
scattering position p 2 is then given by Eq.(fT0|). With a 
branching ratio bfj am ~ 0.15, the scattering process is 
Raman, in which case we assume that the Raman pho- 
ton escapes from the region without further interaction. 
If the scattering process is Rayleigh, then the photon 
propagates along the new direction k G , which is chosen 
from the Rayleigh phase function (e.g. Schmid 1996, Lee 
& Lee 1997a). 

3.3. Ionization Fronts 

The part of the stellar wind from the giant facing the 
hot white dwarf will be ionized by strong far UV radia- 
tion. Taylor & Seaquist (1984) investigated the ioniza- 
tion structures that depend the density and the ionizing 
luminosity. They introduced a a parameter X defined by 



X 



4:TTp 2 m,p 



aL ph { Voo /Mf 



(21) 
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Fig. 2. — Ionization structure in a symbiotic star for a range 
of ionizing luminosities represented by the parameter X defined in 
Eq. (12 1 [ i. The dots at the origin and at the point (80, 0) stand for 
a mass losing giant and a hot white dwarf. The volume of the H I 
region increases as X decreases. The ionization front corresponding 
to X = 0.02 is approximated by a hyperbola given by y 2 = 6.5(ie — 
40) 2 - 400. 

Here, a is the recombination coefficient, L p h is the lumi- 
nosity of ionizing radiation, and a is the photoionization 
cross section. 

Then, the ionization front is described by the equation 



X = f(u,9) 



n(r) s ds 



(22) 



where s is the distance measured from the white dwarf. 
Taylor & Seaquist (1984) provided the closed form of 
f(u, 9) using the density law n(r) oc r~ 2 , which is slightly 
different from the one given by Eq. ([8]) adopted in this 
work. Although the description of Taylor & Seaquist 
(1984) is also approximately valid for the current work, 
we present the detailed ionization structure in a closed 
form in this subsection for the sake of completeness. 

With the density law given by Eq. ((5]) the necessary 
integral to be carried out is 



(pi cos (j) ± \J p 2 — b 2 ) 2 dp 



(23) 



where p and s are related by the impact parameter b as 
inEq. CEU). 

For b = e -1 > 1, an explicit expression of the integral 
I is given by 

Pi cos 2(f> esc (j> tan - 1 y/ (ep) 2 
4e 2 cot 2 6 - 2e 2 + 2 esc 2 A 



I 



: 0) v /eV 



± 2p.i cos (j) ( In 



e 2)3/2 
P 



■ tan 



e(p-l)(l-e 2 ) 
'ep-l 



ep+l 



1/2 
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Fig. 3. — Angle-averaged Monte Carlo line profiles of Raman 
scattered 6825 from Keplerian disks for various Ri/R with Ri 
and R being the inner radius and the outer radius of the disk. A 
stationary scattering H I region is assumed. Main peaks are found 
at ±D , the velocity at the outer rim of the disk, which is set to 
30 km s _1 . 

somewhat larger than the value proposed for V1016 Cyg 
by Jung & Lee (2004). Nevertheless, we perform Monte 
Carlo calculations using this value, because the exact 
value of the half opening angle affects only the overall 
flux of the Raman 6825 feature and the profile shape is 
less sensitive. 



4. RESULTS 



4.1. 



1 



P- 



1 



1 



(24) 



A more detailed calculation is given in the appendix. 

In Fig. [21 we show the ionization fronts for various 
values of X. The horizontal axis coincides with the line 
connecting the two stars. The dot at the origin stand for 
the mass losing giant. We place the white dwarf at the 
point (80, 0), considering the binary separation ~ 80 AU 
proposed by Brocksopp et al. (2004). The ionization 
front corresponding to X = 0.02 is fitted by a hyperbola 
given by y 2 — 6.5(x — 40) 2 — 400. The angle between 
the horizontal axis and an asymptote is 69°. This is 



Pure Keplerian Disk 

In Fig. [31 we show the profiles of Raman scattered 6825 
formed in a stationary H I region relative to a Keplerian 
disk. The H I region is modeled as the hyperboloid given 
in Fig.[2]of which the asymptotic cone has an half opening 
angle o = 69°. In the figure, we fix the velocity v = 
30 km s -1 at the outer edge r — R , and varied the 
inner radius Ri from 0.01i? o to 0.9i? o . In the absence 
of a bright spot, the emission profile is symmetric with 
respect to the line center. 

It is quite clear that the main peaks correspond to the 
velocity v at the outer rim of the disk, which dominantly 
contributes to the whole line flux. The emission region is 
almost ring-like in the case Ri — 0.9R o , where the profile 
has a sharp edge at ±w . With decreasing Ri the profile 
broadens around the sharp edges at ±w . However, when 
Ri < 0.1R o , the profiles are effectively similar, because 
the inner disk is simply physically too small to contribute 
to the overall profile significantly. In view of this, we fix 
the ratio of the inner and outer radii to be Ri/R = 0.1 
from now on. 

In Fig. 21 we show three line profiles for a = — 1, and 
1, where a is defined in Eq. [3] In this figure, Ri — 0.1R o 
is fixed in order to focus on the dependence of profiles 
on the functional form of the emissivity. With negative 
a, the contribution from the inner disk region becomes 
larger, which leads to a broader profile with enhanced 
wings. However, the locations of the main peaks are 
fixed at velocities v — ±w„, yielding qualitatively overall 
similar profiles irrespective of the detailed values of a. 
From this consideration, we fix a — in the remainder 
of this paper. 

4.2. Profile Modulation by the Giant Wind 

In Fig. 03 O VI line photons are incident on the stellar 
wind around the giant that is expanding according to 
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6S15 6820 6825 6830 6835 6840 
Wavelength [A] 

Fig. 4. — Line profiles from a Keplerian disk with the emissivity 
given by e(r) = Ar a , where a = ±1 and a = 0. Solid line shows 
a profile for a = 0, dotted line for a = — 1 and long dashed line 
for a = 1. Far wings, contributed from the inner disk region, get 
stronger as a decreases. 




Fig. 5. — Angle-averaged line profiles of Raman scattered 6825 
from a Keplerian disk that are formed in a spherically expanding 
H I region around the Mira component characterized by the termi- 
nal wind velocity Vca ■ The H I region is partially ionized with the 
ionization front given by the hyperboloid shown in Fig. [2] With 
smaller the double-peak structure is more conspicuous, and it 
is completely erased when «„ > 15 km s — 1 . 

Eq. ©. In this figure, the wind terminal velocity Woo 
has a range from 5 km s _1 to 15 km s -1 . Here, we also 
fix the inner radius Ri = 0.1R o , and other parameters 
to be the same as those considered in Fig. [3] Therefore 
Fig. [3] is obtained when we set Uqo = 0. 

As is seen clearly in the figure, the sharp edges near 
v = ±v are smoothed by the expanding motion of the 
H I region. When Voo exceeds 13 km s _1 the double peak 
structure can hardly be noticed. From this, it is con- 
cluded that the giant stellar wind with terminal velocity 
< 13 km s _1 is consistent with the double peak profiles 
observed in V1016 Cyg and HM Sge. 

However, the profile modulation from other sources 
such as turbulence in the emission region may give rise 
to similar profiles shown in Fig. [5j We do not consider 
this possibility in the current work. 

4.3. VI Resonance Scattering Region 

It is clear that the redward asymmetry exhibited in 
the Raman 6825 features of V1016 Cyg and HM Sge can 
not be accounted for by a pure Keplerian disk emission 
model. It is quite uncertain whether we arc seeing an 
excess in red emission or a deficit in the blue part, or 
possibly both. An excess red emission may correspond 
to a hot spot in an accretion disk which will be discussed 
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Fig. 6. — Best fitting profiles from our Monte Carlo calculations 
and the BOES spectra of V1016 Cyg and HM Sge. Dots show the 
BOES spectra and solid lines represent our Monte Carlo results. 
The fitting parameters for V1016 Cyg are v = 30 km s -1 ,Voo = 
11 km s~i,v c = -10 km s _1 and Av c = 9 km s _1 . For HM Sge 
they are v = 26 km s -1 ,^^ = 10 km s — ,v c = — 7 km s _1 and 
Av c = 6 km s — 1 . 

in the following subsection and in this subsection we con- 
sider the first possibility. 

In order to account for the stronger red part in the 
observed profiles, another component need to be added 
in the previous model considered in Fig.O Therefore, we 
introduce an O VI region between the white dwarf and 
the giant. This additional O VI region is assumed to have 
a very low density so that the O VI emission from this 
region is negligible compared to the Keplerian emission 
disk. Nevertheless, the line center optical depth exceeds 
unity in a velocity space centered at v c < 0. Schmid et 
al. (1999) investigated O VI 1032 profiles and compared 
with those of non-resonant O VI lines to conclude that 
the line transfer effect due to resonance scattering of O VI 
1032 is present in symbiotic stars including V1016 Cyg. 

More specifically, we assume that the scattering optical 
depth by the hot wind component is given by a Gaussian 
function 

t w = exp[(> - v c ) 2 /Av 2 c }. (25) 

This prescription hinders the incidence of those line pho- 
tons with the Doppler factor near v c /c upon the scat- 
tering region. However, the velocity modulation due to 
the giant wind fills the gap significantly, allowing one to 
obtain a double-peak profile with some suppresion in the 
red part. 

Fig. [6] shows the best fitting profiles to out BOES 
spectra of V1016 Cyg and HM Sge. For V1016 Cyg, 
the velocity v at the outer rim of the accretion disk 
is v a = 30 km s _1 and the terminal velocity of the gi- 
ant wind Voo = 11 km s . The corresponding val- 



ues for HM Sge are v = 26 km s 



10 km s" 



The parameters adopted for the hot wind component in 
V1016 Cyg are v c = —10 km s _1 and Av c — 9 km s _1 
and v c = —7 km s _1 and Av c — 6 km s _1 for HM Sge. 

In Fig. [7] we present the line profiles for various veloc- 
ities v at the disk outer edge ranging from 24 km s _1 
to 33 km s _1 with Voo,v c and Av c fixed to the best fit 
values of Fig. ([6|). The middle panel shows the best fit 
line profiles in Fig. [SI whereas the top and bottom pan- 
els show line profiles with faster and slower disks, re- 
spectively. Considering the much poorer fits shown in 
the top and bottom panels, we may conclude that the 
outer rim velocities are constrained to be in the range 
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FlG. 7. — Line profiles from our Monte Carlo calculations for 
various velocities vq at the disk outer rim. Dots show the BOES 
spectra and lines represent our Monte Carlo results. Solid lines 
show our best fit profiles shown in Fig. 6. 
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Fig. 8. — Line profiles from our Monte Carlo calculations for 
various terminal giant wind velocities Voo with v , v c and Av c fixed 
to the best fit values of Fig. JB}. Dots show the BOES spectra and 
lines represent our Monte Carlo results. The best fit profiles in 
Fig. 6 are shown by solid lines. Larger values of Voo tend to yield 
better fits at far wing parts but erase the double peak structure in 
the main parts. 

v = 30±3 km s- 1 for V1016 Cyg and v Q = 26±2 km s" 1 
for HM Sge. Adopting a typical white dwarf mass of 
Mwd = 0.7 M©, the main O VI emission region resides 
from the hot white dwarf component at ~ 0.7 AU for 
V1016 Cyg and - 0.9 AU for HM Sge. 

Fig.[8]shows the Monte Carlo line profiles computed for 
various terminal giant wind velocities Voo from 8 km s _1 
to 13 km s _1 also with v Q , v c and Av c fixed to the best fit 
values introduced in Fig. ([B]). Similarly as in Fig. the 
middle panels show our best fit profiles in Fig. [6] for com- 
parison. As is discussed in the previous section, larger 
giant wind velocities tend to erase the double peak struc- 
ture, whereas overall sharper profiles result from smaller 
wind velocities yielding poor fits at the outer parts of 
the observed profiles. It is also interesting to note that 
the wind terminal velocity Voo ~ 10 km s _1 for both 
V1016 Cyg and HM Sge is quite similar to the escape 
velocity of a typical giant star. 

4.4. Hot Spots 

As is well known from spectroscopic studies of cata- 
clysmic variables, an accretion disk may have a bright 
spot which is formed by the impact of the accretion 
stream near the outer edge of the accretion disk (e.g. 
Warner 1995). The existence of a similar feature in sym- 
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FlG. 9. — Monte Carlo line profiles from a Keplerian disk with 
a hot spot. The solid lines show our Monte Carlo data and the 
dots show the BOES data. The dotted lines show the best fit 
profiles given in Fig. [6] The fitting parameters are v = 33 km s — 1 
and Voo = 9 km s — for V1016 Cyg and v a = 27 km s _1 and 
Voo = 9 km s _1 for HM Sge. 

biotic systems is not clear. In this subsection, we inves- 
tigate the effect of the introduction of a hot spot in a 
Keplerian disk on the line profile. 

In our Monte Carlo code, photons generated in the 
specified region satisfying r > 0.8r o and 10+0. 77r| < 0.17T 
is assigned a statistical weight of 4.8, so that this region is 
4.8 times brighter than other part of the emission region. 
This mimics the bright spots well known in cataclysmic 
variables, for which the accretion process takes place via 
Roche lobe overflow. Here, we do not consider the hot 
wind component that provided central absorption part 
in Fig. El 

Fig. shows our profiles. The solid lines show the 
profiles from a Keplerian disk with a hot spot, whereas 
the dots show the BOES data. For the sake of compar- 
isons, we also show the best fit profiles given in Fig. [S] 
by dotted lines. The fitting parameters for V1016 Cyg 
are v = 33 km s _1 and Voo = 9 km s _1 . For HM Sge, 
they are v — 27 km s^ 1 and v^o — 9 km s _1 . It should 
be noted that the disk outer rims move slightly faster in 
this model than the ones considered in the case of the 
best fit model. 

The profiles provide good fits to the outer parts A < 
6824 A and A > 6830 A, whereas poor fits are obtained 
near center. The introduction of a hot spot enhances 
the red peak. However, the conspicuous central dips in 
the observed data are not obtained in pure Keplerian 
disk models with v > 25 km s _1 for Voo > 8 km s _1 . 
Lowering makes the fit worse by narrowing the widths 
of the two peaks. 

5. SUMMARY AND DISCUSSION 

In this work, we present our profile analysis of Raman 
scattered O VI 6825 in the symbiotic stars V1016 Cyg 
and HM Sge, assuming that O VI line photons are gener- 
ated from a thin Keplerian disk with a scattering region 
between the white dwarf and the Mira giant. The neutral 
scattering region is assumed to be a spherical slow stellar 
wind with terminal wind velocity with an ionization 
front approximated by a hyperboloid. The best fitting 
parameters are v = 30 km s _1 , = 11 km s _1 and 
v w = —10 km s _1 for V1016 Cyg and v = 26 km s _1 , 
Voo — 10 km s^ 1 and v w = —7 km s _1 for HM Sge. 

The overall fitting quality is quite insensitive to our 
choice of the inner radius of the O VI emission region 
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and the functional dependence of emissivity on radius. 
Furthermore, the profile shape is little affected by the 
exact shape of the ionization structure, which is deter- 
mined by the mass loss rate of the Mira component and 
the ionizing luminosity from the white dwarf component. 
However, the ionization structure mainly affects the total 
flux of Raman scattered 6825. In this point, an indepen- 
dent study should be added to current work in order to 
adequately probe the ionization structure and mass loss 
and transfer processes in symbiotic stars (e.g. Jung & 
Lee 2004). 

Although the introduction of a hot spot in our Keple- 
rian disk model failed to improve the profile fitting qual- 
ity in the case of V1016 Cyg and HM Sge, this should 
not mean to exclude the possibility of the asymmetri- 
cal emission in symbiotic stars. As Harries & Howarth 
(1996) reported, triple peak structures in Raman scat- 
tered 6825 in symbiotics are more prevalent than double 
peak profiles in symbiotic stars. It is quite probable that 
more complicated profiles may be fitted using more in- 
gredients in the emission and/or scattering regions. 

The resonance scattering O VI component considered 
in section 4.3 may be envisioned as a part of ionized giant 
stellar wind, in the vicinity of the inner Lagrange point. 
In order to scatter a significant fraction of slightly blue 
O VI photons, the region should move toward the main 
O VI line emission region around the white dwarf. We 
may note that the bestfit speed v c in Fig. [5] is slightly 
less than the giant wind terminal velocity. Because v c 
is the velocity component along the direction connecting 
the giant and the white dwarf, we may interpret that 
the resonantly scattering O VI region passes through the 
inner Lagrange point making an angle 

6, ~ sin" 1 0.8 = 50° (26) 

under the assumption that the speed of the region is co- 
incident with the wind terminal velocity Woo- With the 
wide binary orbit of 50 — 80 AU the orbital speed of the 
giant component in both V1016 Cyg and HM Sge is com- 
parable to v c . Although we only consider the spherical 
wind around the Mira type giant component in the cur- 
rent work, it is also a possibility that the wind around the 
giant component also possesses an azimuthal component, 
in which case this inference of 0i should be corrected ap- 
propriately. Refined hydrodynamical studies may shed 
more light on the mass transfer processes in symbiotic 
stars. 

The red wing excess in V1016 Cyg that is not well fit- 
ted by our model may indicate the existence of another 
emission component and/or a neutral scattering region. 
Spectropolarimetric studies by Harries & Howarth (1996) 
show that the far red wing part is polarized in the direc- 
tion perpendicular to that for the remainder part. These 
facts are consistent when the additional emission and/or 
neutral component is receding with a velocity of several 
tens of km s _1 perpendicular to the binary orbital plane. 

It is unclear that this component is associated with the 
bipolar structure of these objects. Solf (1983, 1984) de- 
tected a bipolar outflow in [N II] lines in HM Sge. The 



velocity scale of several tens of km s _ 1 is much smaller 
than the usual velocity scale associated with the fast hot 
wind that may emanate from a deep inner region of an 
accretion disk in the vicinity of the hot star. Schmid et al. 
(1999) reported the existence of broad wing components 
in O VI 1032, 1038 observed with ORFEUS. Even though 
the broad wings are consistent with the electron scatter- 
ing origin, a fast bipolar wind around the hot component 
can also be a candidate. It is an interesting possibility 
that the component responsible for resonance scattering 
of O VI is a clumpy one driven by the fast hot wind 
moving in the polar directions. More interesting results 
arc expected when photoionization calculations and hy- 
drodynamical studies can be combined with spectropo- 
larimetric polarimetric observations of these symbiotic 
objects. 

It is also interesting that these two symbiotic stars 
exhibit He II Raman scattered features blueward of H 
Balmer lines. In particular, Raman scattered He II 4850 
features of V1016 Cyg and HM Sge exhibit no multiple 
peak structures (e.g. Jung & Lee 2004, Birriel 2004). 
Because the cross section for He II Raman scattering is 
larger than that for O VI Raman scattering by two orders 
of magnitude, we may expect that the scattering region is 
much more extended than the Raman O VI counterpart. 

The profiles of Raman O VI 6825 in V1016 Cyg and 
HM Sge are totally different from other high ionization 
lines including He II. For example, He II 4686 emission 
lines of V1016 Cyg and HM Sge are single-peaked with 
slight blue asymmetry. Robinson et al. (1994) consid- 
ered line profiles from accretion disks in symbiotic stars, 
and it appears to be difficult to infer the existence of 
an accretion disk in symbiotic stars from these emission 
lines. It may be that the He II emission region is more 
extended than the O VI emission region. It is beyond 
the scope of this work to model other emission lines and 
more sophisticated models including photoionization and 
hydrodynamics are required. 

The disparity of the profiles in Raman scattered O VI 
6825 and Raman scattered O VI 7088 can be an impor- 
tant clue to the structure of the accretion disk in sym- 
biotic stars. A number of researchers pointed out the 
relative weakness in the blue part of Raman O VI 7088 
(e.g. Schmid et al. 1999, Harries & Howarth 1996). The 
profile difference in resonance doublets of Si/ 2 — -Pi/2,3/2 
was reported in symbiotic stars and young planetary neb- 
ulae (e.g. Feibelman 1983). 
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APPENDIX 
THE IONIZATION FRONTS 

We describe in detail the photoionization calculation for the ionization fronts, which was briefly discussed in section 

2. 

The distance from the giant is measured in units of the radius of the giant, by p = r/R*. Let p t be the separation 
of the giant and the white dwarf. Along a path from the white dwarf making an angle <f> with the line connecting the 
giant, the physical distance is measured by the parameter s satisfying 



s = 



_ J p.i cos (f) — \J 'p 2 — b 2 for s < pi cos < 
Pi cos <j> + \J p 2 — b 2 for s > pi cos < 



(Al) 



where b — pi sin <j> is the impact parameter from the giant to the path of the ionization radiation from the white dwarf. 
With the density law 

"<"» = *rry < A2 » 

where the distance r from the giant center is measured in units of i?* by the parameter p — rjR^. 

The recombination rate along the path of ionization radiation originating from the white dwarf is proportional to 
the density squared along the path, given by 



Jn(p) 



2 s 2 ds. 



(A3) 



In terms of p, this integral can be re-written as 



(ft COS0± y/p 2 - b 2 f 
p{p-\) 2 ^ 2 ~b^ 

= pj cos 2 <t>I 1 +I 2 ±2p i I 3 , 



dp 



(A4) 



We define the three integrals as 



h 

h 
h 



I 



dp 



I 



p( p-i) 2 ^p^¥ 
^P^Vdp 
p(p-i) 2 

dp 



p(p-l) 2 - 



The integral -Z3 is elementary and is given by 



Considering the substitution 

the integrals under consideration become 

1 



p-l p-l 
p = b cosh 



(A5) 

(A6) 
(A7) 



h= j d9b 2 (cosh 2 8- 1) 



+ 



1 



b cosh 8 b cosh 8—1 (frcosh# — l) 2 
1 1 



1 



b cosh 8 b cosh 8—1 (b cosh 8 — 1) 
In order to obtain the closed form expressions of the integrals Ii and I 2 , we define the two integrals K and L as 

d8 



(A8) 



L- 



cosh 8 — e 
d8 

(cosh8-e) 2 ' 



(A9) 



with € — 1/6. 

For e < 1, making use of the half angle formula for cosh# and the substitution 



tanh - = 



1 - e 
TT~e 



u, 



(A10) 
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we obtain 



K = 



: tan 



1 



1 



■ tanh ■ 



On the other hand, for e > 1, then the substitution we consider is 

, 9 
tanh — = 
2 

which leads us to 



e- 1 



K 



dv 



Ve^T J v 2 - I 
We also consider the following relation 

d ( sinh 9 



coth 



— tanh — . 
1 2 



from which we obtain 



d9 \ cosh 9 — e 



L 



1 



1 



rK 



cosh6> — e (cosh6* — e) 2 ' 
1 sinh0 



1 — e 2 cosh 9 — e 



The integrals I\ , Ii in Eq. (|A8j) can be finally expressed as 



h 



1 



• tan sinh 9 



K 

T 



L 



1-2 = —b tan sinh 9 



(Ir 



b b 2 

Combining these results, for e < 1, an explicit expression of the integral / is given by 



I = pi cos2(/>csc</Han 1 sinh 9 
4e 2 cot 2 6 - 2e 2 + 2 esc 2 M 



esc 



sinh f 



e(l 

± 2pi cos 4> ( In 



e 2)3/2 

p 



■ tan 



1-e 2 
-l 



cosh 9 — e 



1 



1 



— tanh — 
€ 2 



1 



Also, for e > 1, we have 



1 



I — Pi cos2</>csc0tan sinh6> 
4e 2 cot 2 d> — 2e 2 + 2 esc 2 6 



esc 



sinh 9 



e(e 2 - 
— 2pi cos 4> ( In 



!)3/2 

p 

o-l 



1 - e 
coth -1 



cosh 9 



1 



tanh — 
1 2 



1 



(All) 

(A12) 
(A13) 

(A14) 
(A15) 



(A16) 



(A17) 



(A18) 



Noting Eq. (|A7p . we may immediately see that Eq. (|A17[) reduces to Eq. (|2"4"]) . In Fig. [21 we show ionization fronts 
for a range of ionizing luminosities characterized by the parameter X. The cool giant is located at the origin and the 
white dwarf is at (pi — 80, 0). We also approximate the ionization corresponding to X = 0.02 by a hyperbola given by 



V 



6.5(a; - 40) 2 - 400. 
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